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, The zero-temperature dynamical spin-spin correlation functions are calculated for the spin-| 

two-leg antiferromagnetic Heisenberg ladder in a magnetic field above the lower critical field H c \. 
0^ ' The dynamical structure factors are calculated which exhibit both massless modes and massive ex- 

citations. These modes appear in different sectors characterized by the parity in the rung direction 
5— i ■ and by the momentum in the direction of the chains. The structure factors have power-law singular- 

ities at the lower edges of their support. The results are also applicable to spin-1 Heisenberg chain. 

■ The implications are briefly discussed for the various correlation functions and the 7r-resonance in 
the SO(5) symmetric ladder model. 
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£j ' I. INTRODUCTION 
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Systems locating in between one dimension and two dimensions, ladders, have attracted great attention both 
theoretically and experimentally.!!] This line of research has been pursued both with the hope of finding some clues 
to the understanding of the high-temperature superconductivity, and because of the experimental discoveries of 
new ladder compounds, such as a cuprate spin laddeia SrCu2C>3 and a superconducting laddeiu Sr a Cai4_ x Cu2404i. 
Another spin ladder system of current interest is the organic material Cu2(C5Hi2N2)2CLi whose spin-gap behavior ■ 
observed by the measurements of the spin susceptibility, the magnetization curve, and the NMR relaxation rate! 

i ^ i ■ The nonzero magnetization is observed once the external magnetic field H exceeds the lower critical field H c x which is 
equal to the spin gap (in units where the product of the Bohr magneton and the (/-factor is set to be unity). This is a 

C*~) . quantum phase transition driven by the external magnetic field from the gapped spin liquid to the gapless Tomonaga- 
Luttinger liquid state. These experiments motivated further theoretical w-acks on the various properties of the. gapless 
regime: the magnetization process-pLl the spin-spin correlation functions, mil and the spin-Peierls instability-ErE-il 
In this paper we extend the theory of Shelton et al&3 to the two-leg Heisenberg ladder in the gapless regime, 

\ (H > H c \) and calculate the dynamical spin-spin correlation functions and structure factors in the ground state.EJ 

We obtain the dynamical structure factors containing massive spin excitations as well as massless excitations. The 
latter cont ribution is characteristic of the Tomonaga-Luttinger liquid and is commonly found in the S — | Heisenberg 

■ chainiSB This component may be interpreted as coming from the bose-condensate of S z — 1 magnons .E3 On the 
other hand, the massive excitations we are interested in have S z — magnons in their origin. Below H c \ the massive 
magnons give rise to (5-function peaks in the dynamical structure factor. In the gapless regime (H > H c i) the 5- 
function peaks of the massive excitations turn into power-law singularities, because the massless excitations introduce 

, algebraically decaying prefactors to the exponentially decaying correlation functions. An interesting feature of the 
antiferromagnetically coupled Heisenberg ladder is that these massless and massive modes appear in different sectors 
of the structure factors, characterized by the momentum q y in the direction of the rungs (q y can take only and it) 
O ■ and by the momentum q in the direction of the chains. 

Oup-study of the massive excitations in the gapless regime is also motivated to better .understand the 7r-resonance 
•*h , modet3 in the SO(5) symmetric ladder model introduced recently by Scalapino et alxJi In this model there is a 
quantum phase transition driven by the chemical potential from the spin-gap Mott insulator to the d-wave-like 
superconducting phase. It was shown in the strong-coupling limit that the spin-gap magnon mode of the Mott 
insulator evolves continuously into the 7r-resonance mode of the superconducting phase.E3 A complementary weak- 
coupling acproachEScJ was taken to study a general two-leg ladder model of weakly interacting electrons. In particular, 
Lin et altB showed that the model is renormalized to a fixed point where a global SO (8) symmetry is realized that 
contains the SO (5). They also obtained the ground-state phase diagram and low-energy excitation spectra at half 
filling which are qualitatively in agreement with the strong-coupling picture of the SO (5) symmetric ladder model. Yet 
the energy and the spectral weight of the 7r-resonance mode in the dynamical spin structure factor are not completely 
understood even for the ideally constructed SO (5) symmetric ladder. An extensive numerical exact-diagonalization 
study is performed to clarify this and related issuesJ^J but clearly it is desirable to develop an analytic theory to 
discuss the power- law singularity which is expected to appear in the structure factor O The field-driven quantum 
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phase transition in the Heisenberg ladder we study in this paper is analogous to the quantum phase transition in the 
SO (5) symmetric ladder. A S z = 1 boson in the Heisenberg ladder corresponds to a hole pair in the latter model, 
and a S z — boson to a magnon or the 7r-resonance mode. Therefore our Heisenberg spin ladder may be viewed 
as a toy model for the SO(5) symmetric ladder. Using the analogy, we can deduce the correlation functions for the 
d-wave superconductivity, the charge density wave, and the spin correlations in the SO (5) ladder from the spin-spin 
correlators of the Heisenberg ladder. 

This paper is organized as follows. In Sec. II we introduce the model and analyze it in the limit of strong interchain 
coupling, where it is easier to get intuitive pictures of the physics. Detailed calculation of the dynamical spin-spin 
correlation functions and the structure factors are presented in Sec. III. Here we take the opposite limit of the weak 
interchain coupling, and employ the bosonization method. In Sec. IV we briefly discuss implications to the 7r-resonance 
in the SO (5) symmetric ladder model and summarize the results. 



II. STRONG-COUPLING APPROACH 



The Hamiltonian for the two-leg Heisenberg ladder we study in this paper is 

oo oo oo 

h = j Y s ^- i ■ s ^- i + i + J± Y ■ s 2 ,t -h s* >u (i) 

/i=l,2 i— — oo i— — oo fj.= l 7 2 i— — oo 

where S^i is spin-i operator, and the intrachain coupling J is positive (antiferromagnetic) . The interchain coupling 
Jx is also assumed to be positive, unless otherwise noted. 

In this section we briefly discuss static correlations above the lower critical field H c \ in the strong-coupling limit, 
Jx » J. When J/J± = 0, the ladder is decomposed into independent rungs, each rung consisting of two spins 
Sij and S , 2,i- There are four eigenstates in each rung, the singlet \Si) = (|tl) — |lt))/V2) an d the triplet states 
K+) = | TT), \T h0 ) = (I Tl) + | IT»/V2, and |T<_) - 1 Their energies are E(Si) = -3J ± /4, E(T i>+ ) = J ± /4 - H, 
E(Ti t o) = Jj_/4, and E{Ti.-) = Jj/4 + H, respectively. It is then natural to introduce boson operators s\, ij+> 



t| , and t\ _, which create the singlet state and the triplet states at the ith rung.E3Sla They obey the constraint, 
s l s i + 4 + 4 o^,o + 4 -U,- = 1- With these hard-core bosons we represent the spin operators 

(2a) 
(2b) 
(2c) 

(2d) 

(2e) 

(2f) 

The suffixes and it refer to the momentum q y . Note that Sq i and do not contain the t^o boson. 

In the presence of a high magnetic field H > Jj_, we may ignore the i^o and t^- bosons as they have higher energies 
than the boson. In lowest order in J the effective Hamiltonian thus becomes 
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with the simplified constraint s|s, + 4 = 1- As pointed out in Ref. 24, this effective Hamiltonian may be written 
as 

H cff = JJ2 (sfSf +1 + SfSf +1 + \siSt + ^j + (jx + |-#)£^ + const < W 
where is a spin-i operator 
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S+ = t[ + s h (5b) 
Sr = s \t,. + . (5c) 

We notice that Eq. (|j) is just the Hamiltonian of the spin-i XXZ chain with the external magnetic field H = 
J± + j—H, whose properties are well understood. When \H\ < | J, the XXZ chain is not fully polarized. El This means 
that the ladder has unsaturated mag netization (5^) = = M (0 < M < 1) when J± - J < H < J± + 2 J. 

We have thus obtained the lower and upper critical fields H c \ = J± — J and H C 2 = J± + 2J in lowest qwder in J/Jl.l£| 
The magnetization curve for M < 1 can also be obtained from the known result for the XXZ chaincJ 



+ -^p!£M. (6) 



A similar result holds for 1 — M <C 1. The square-root behavior is a well-known universal behavior .lZTlj'EI When 
< M < 1 the ground state is a superfluid of the bosons. Within our approximation the following equal-time 
spin correlation functions are readily obtained from those of the XXZ chainEd 

(Sln+A) = ((S* n+ i + V2)(5f + 1/2)) = M 2 - * + ^ cos(27rMri), (7a) 

(5*„ +1 5* x ) = 2<5* +1 5*) = c 3 - c 4 |^ cos(27rMn), (7b) 

where c's are numerical constants. For M<1 the parameter r\ is given b 

T) = 2— -M + 0(M 2 ). (8) 

In general < 2 for < M < 1, and 77 — > 2 when M — + 0, l.H Thus, the leading term in (<Sq n+1 /Sg 1), besides 
the constant M 2 , is n - '' cos(27rMn), which reflects the fact that the hard-core bosons t i + have a tendency to form 
a density wave with a period equal to 1/M. On the other hand, the correlation functions (S*^S*j) and (Sq^Sqj) 
involve the massive tifi and bosons, and therefore decay exponentially for |z — j| 3> 1. We conclude that the 
correlators (SqSq) and (S^S^) show quasi-long-range order while (S*S*) and (Sjf fiff) are short-ranged. These 
correlation functions are discussed in more detail in the next section. 

Before closing this section, we make a few comments. First, Eqs. ([7a|) and (|7b| ) have the same form as the pepeiators 
(S^ +1 Si) and {S* +1 Sf) of the 5=1 Heisenberg chain in a magnetic field larger than the Haldane gap.E3'EjH An 
important difference is that rj > 2 in the 5=1 chain.E3 We will ccpac-back to this point in the next section. Second, 
from the knowledge of the exact propagator of hard-core bosons P 4 " 35 ! we may expect that the correlation function 

(5J n+ i5^ L ) should have terms proportional to cos(27r/Afn) with I = 0, 1, 2, Third, as noted in Introduction and 

will be discussed in Sec. IV, the boson and t^o boson correspond to a pair of holes sitting on a rung and to a 
magnon in the SO(5) symmetric ladder, respectively. Thus, we may compare the hole-pair correlation and the spin 
correlation functions with (S* iS% *) and (S* •) of our model. 



III. WEAK-COUPLING APPROACH 



In this section we calculate the dynamical spin-spin correlation functions in the weak-coupling limit., Jk. <gj_ J. 
We use the Abelian bosonization method which has been successfully applied to the Heisenberg ladder .lSEZHIj'EI In 
this approach we first bosonize two independent spin-i Heisenberg chains and then treat the interchain coupling J± 
perturbatively. It is a relevant perturbation and is renormalized to a strong-coupling regime, generating a mass gap 
in the excitation spectrum. This behavior does not depend on the sign of Jj_, and therefore the model describes both 
the antiferromagnetic Heisenberg ladder for Jj_ p*. and the 5=1 Heisenberg chain for J± < 0.E3 The spin gap in 
the latter case is nothing but the Haldane gap.E3 Since the gapless phase in a magnetic field is a single phase for 
< Jj_ < 00, we expect the correlation functions in the weak-coupling limit to have the same structure as in the 
strong-coupling limit. ■— ■ 

We follow the formulation of Shelton et aZ.,t3 which we explain below to establish the notation. We begin with the 
bosonized Hamiltonian of the ladder in the continuum limit. It consists of three parts: TL = TL+ + 7Y_ + Ti±, where 
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H+= I dxl- 
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(9) 
(10) 

(11) 



Here w is the spin- wave velocity, do is a lattice constant, and m = J±\ 2 /2ir with A being a numerical constant. 
The bosonic fields <j>± and 9± obey the commutation relations [cf) + (x),9 + (y)] — [(f)- (x), 9- (y)] = iQ(y — x) and 
[<f) + {x) 1 4>-{y)} — [cj>+(x),8-(y)] = [<f)-(x), 9+(y)] = [9+(x),8-(y)] — 0, where 8(x) is the step function. It is important 
to note that H appears in TL + only. Thus, the external magnetic field changes the dynamics of the fields (/>+ and 9+, 
while the other fields cj>- and 9- are not directly influenced by the uniform field H. The excitations involving these 
latter fields remain gapful even above the lower critical field H c \. The spin operators with q y — and tt, defined in 
Eqs. (pa|)-(p^), are written in terms of the bosonic fields as 



S*(x) 
S+(x) 



Qq d<j) 
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2X 



t/tt dx 
Oq dip- 



(-l) l/fl »-sin(^ + )cos(^_), 
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dx 



(_l)x/ao — C0S(^+) Sin(^- 
7T 



(12a) 
(12b) 



S+(x) = 



(— l) x / a °\cos(</n9-) + cos(y/n9-) sin( v / 7r0+) cos(\/7r</>-) + % sin(V7r#_) cos(\/7r^ + ) sin(v / 7r</)_) , 

(12c) 

i(-l) x / a °Asin(v / 7r^-) + cos(Vtt^-) 008(1/7?^+) sin(A/7T</)_) + zsin(-v/7T0-) sin(v / 7T0+) cos(-v/7r0-) 



(12d) 



for x 



iao (a = x,y,z). Equations (12aJ)-(12d) directly follow from Eq. (30) and Appendix A 



where 5° (x) 
of Ref . [lj. 

Shelton et al. showed that TL+ and 7i_ can be greatly simplified by fermionization. The Hamiltonian 7i_ then 
becomes that of free massive Majorana fermions £3(2;) and p(x) having the mass gap m and —3m, respectively. This 
is equivalent to the two-dimensional Ising model above or helow the critical temperature. This observation allowed 
them to obtain the dynamical spin-spin correlation functionsLj from the known results of the Ising model.Ell Physically 
the £3 fermion describes the S z — magnon excitation, whereas the p corresponds to a singlet excitation with much 
higher energy. 

On the other hand, the <f>+ and 9 + fields represent the S z = ±1 magnon excitations. As the S z = 1 bosons 
condense above H c i, these bosonic fields have massless excitations. We will first integrate out the massive Majorana 
fermions and concentrate on the massless bosonic fields. To proceed, here we introduce an approximation for dealing 
with Ti.±. This interaction Hamiltonian has three components. The first component involving only </>_ and i.e., 
cos(V47T0_) cos(V4tt0-) and (d(f>- j dx) 2 , has two major effects on the dynamics of H-. Ooe effect is to renormalizc 
the bare mass m to mln(A/m), where A is a high-energy cutoff, as noted by Shelton et a/.l!3 This can be absorbed by 
redefining the mass. The other effect is a strong twpTparticle collision described by a S matrix having a superuniversal 
form, as recently discussed by Damle and Sachdev.cil Since we are only concerned with processes in which at most one 
S z = magnon is created, this strong scattering effect may be irrelevant for our discussion of the dynamical correlations 
at zero temperature. The second component is a coupling term, cos(V4~7r6'_) cos(-\/4~7r<^+)- When integrating out the 
field perturbatively, we find that the leading term (cos(\/47r#_))_ cos(V4~7r0 + ) cx ln(A/m) cos(v47r</ > +) gives the 
renormalization of the mass of the <fi + and 9+ fields, m — > mln(A/m), as expected from the SO(3) symmetry, where 
the average is taken in the ground state of 7i_. This is again taken care of by redefining the mass. The higher-order 
terms will generate, through gradient expansions, both irrelevant terms like cos(-\/4~7rZ0+) with I > 1, which we can 
safely ignore, and a marginal operator (d x <f)+) 2 , which should be kept. The third component is the term (d(j)+/dx) 2 
already present in Ti.± and will be kept in the following calculation. Hence we reduce Tl± to the form 



J±a<3 

47T 



dx 



dx 



(13) 



where the coupling constant J±ao may be modified from its bare value because of the emergence of the operator 
(d x (j> + ) 2 in the higher-order terms discussed above. We believe, however, that the sign of the coupling constant is not 
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changed by the renormalization.B Having made this approximation, we now integrate out the <p- and fields to get 
the spin-spin correlation functions. Within our approximation the fields 0_ and are independent of <p + and #_|_, 
and therefore the correlation functions of 0_ and 6- are independent of H. We use Eq. (33) of Ref. [l2] to represent 
cos(y / 7r0_), sin(y / 7r</)_), cos(y / 7r0_), and sm(y^F#_) in terms of the order and disorder parameters of an Ising model. 
We then use Eqs. (38) and (39) of Ref. O to obtain their correlation functions. The correlation functions (S^Sq 1 ) 
involve cos(y / 7r0_) cos(y / 7r#_) and sin( v / 7r0_) sin(y / 7r0_). They are equivalent to free massive Majorana fermions, £3, 
whose correlators are easily obtained. Finally, the correlator (d x (j)-(x)d x <j)-(0)) decays exponentially [oc e ~ 2mr / v with 
r = (x 2 +u 2 r 2 ) 1 / 2 ] and is ignored. Hence, we arrive at the following expression of the dynamical spin-spin correlation 
functions: 



(^(x,r)^(0,0)) = ^(^<Mx,t)^<M0,0)) + , 

7T 

flW 2 A 2 

(S*(x,t)S^0,0)) = (-l) x/ao — ^A'o(mr/ U )(cos[V^0 + (x,r)]cos[^ + (O,O)]> + , 

V 7T / 7T 

(5+(aj,r)S-(0,0)) = {-If'** (jf) A 2 ^W»+M e -(^ + (o,o)^ , 

(S+(x,t)S q -(0,0)) - _^/ e V5f» + (x,r) e -V5r* + (»,r) ei Vir* + (0 J 0) e -V5f» + (0,0)\ a+iro ( TOr/v ) 
Z7T V \ I + 

a ° / e iV^e + ( a; ,r) e j v ^F0 + ( a; ,r) e - lv ^F0 + (o,o) e -i v ^ : e + (o,o)\ d_K (mr/v) 
2ir 2 v \ / + 



(14a) 
(14b) 

(14c) 



+— -j-iforor/i) 
2ir z v 



e iV^»+ (x,t ) e -iv^F0+ (:t,t) e -iV?</>+ (0,0) e -iV^0+ (0,0) 



e i v /7Fe + (a;,r) e iV7F0 + (a;,T) e i v ^r0 + (O,O) e - 



-iV5Fe + (o,o)^ 



(14d) 



where r is the imaginary time, vl i is a nume rical constant, 9± = d T ± ivc^, and ( ) + is the average with respect to the 
Hamiltonian TL + + H±. In Eqs. ( |l4a|) -(14d) we have ignored the terms decaying much faster than e - mr / v _ Therefore 
we discarded the contribution from processes involving more than one massive magnons. These equations are valid 
for J± > 0. When J± < 0, on the other hand, the strongly renormalized interchain coupling combines the spins 
and $2,1 into a single spin, and the ladder behaves as a S = 1 Heisenberg chain. As explained in Ref. [l^, when taking 
average over the massive Majorana fermions £3 and p, we only need to exchange the order and disorder parameter of 
the Ising model. Using the correlators of 0_ and 6- for J± < in Ref. [l|, we find 



(S^(x, t)S$(0, 0))\ j±<0 = (S * r)S o z (0, 0)) | J±>Q + (Sfa, r)5*(0, 0)) | J±>Q , 
(S+(x, t)S^ (0, 0))| Jx<0 = (S+(x, t)Su (0, 0)) I Jx>Q + (S+(x, r)S~ (0, 0))| J±>Q , 



(15a) 
(15b) 



where the correlators in the right-hand side are those in Eqs. (14a) (14d). When J± < 0, the q y — ir correlators decay 
much faster than the q y = correlator and are thus negligible. Hence, the dynamical spin-spin correlation functions of 
the S = 1 Haldane chain are linear combinations of those q y = and ir correlators of the Heisenberg ladder ( J± > 0). 

Now our task is to calculate the correlators of (j>+ and 8 + in the presence of the magnetic field. The Hamiltonian 
7i+ is identical to the one used to study the commensurate-in comm e nsur ate transition in classical two-dimens; 
systemsJl3~E3 In fact some of the correlation functions in Eqs. (14a)-(14d) have been discussed in this context] 
In particular, the leading term of the correlation function corresponding to (e Jv/ ^+( :E '' r ' 'e^v^+C 3 * )) is obtained by 
Schulzu£l including its universal exponent in the limit M — » 0. These results are used to calculate the spin-spin 
correlation functions by Chitra and GiamarchiJ3 who unfortunately seem to have overlooked some terms including the 
leading term [ex cos(27rMa;)] in (S^(x)S^(O )). We thin k therefore that it is still worthwhile to describe the calculation 
of the correlations of cj>+ and 8 + in Eqs. (14a)-(14d) in some detail, despite the fact that the Hamiltonian H+ has 



been analyzed in many literatures. To our knowledge, our results concerning the massive modes are new. 
Following Ref. O, we fermionize H+: 



TL+ = I dx 



dk 



vk 



(4 



R. k CR,k ~ C L ,k c L,k 



(4 



mi I r), , i'L,k ~ c\ k C B .,k) - H (4 fe Cfl,fc + c[ k C L ,k 



(16) 



where Vl is the right-going (left-going) complex fermion field, and ipRtL)(x) — § {dk/ \f2~Ti)e lkx c R ( L ^ k . The 

fermion fields are related to the bosons by the standard relations: 
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2y / 7T dx 

1 d 



dx 



tj> + (x) - 



f(*)], 

9 + (x)]. 



(17a) 
(17b) 



It is important to note that the normal ordering in the above equations is defined with respect to the ground state of 
H = 0. The fermionized Hamiltonian Ti + is easily diagonalized: 



n+ = 



dk 



\J v 2 k 2 + m 2 — H) alak — (v v 2 k 2 + m 2 + 



> 2 ->-H) 4a fe 



where 



cos(<^fc/2) -isin(^k/2) 
-isin(^ fe /2) cos(v? fe /2) 



(18) 



(19) 



with taxi^fe = m/vk. The magnetic field plays a role of the chemical potential to the fermions. We are concerned 
with the case where H is slightly above the lower critical field H c \(= m) such that < M <C 1. The ground state is 
obtained by filling the upper band (a k ) up to the Fermi points (|fc| < kp) and the lower band (dk) completely; see 
Fig. [l|. The Fermi wavenumber kp is related to the magnetization: fcj? = irM/ao. This follows from 



M 
«o 



dx(d x (j>i 



dx(:ip t R <il) R + ipl'il) L :) + 



1 

27 



dk(a\a k ) 



+) 



(20) 



where L is the length of the ladder and (a k ak)+ = Q(kp~ \k\). In calculating long-distance correlations, we can safely 
ignore the lower band and keep only the low-energy excitations around the Fermi points in the upper band. In the 
fermion representation the interaction term Ti.± reads 



J±a 



dx (f/^VO 



(21) 



where ip(x) = J (dk j \Z2n)e lkx a k and we have dropped the contribution from the lower band. The total Hamiltonian 
for the fermions in the upper band, Ti a = 7i + + Ti± , consists of the kinetic energy, Eq. ( |l8| ) , and the short-range 
scattering term, Eq. (pT|). The coupling constant of the latter term is proportional to J± in lowest order. Thus. 
the interaction is repulsive for the antiferromagnetically coupled ladder, while it is attractive for the S = 1 chain.E 2 ] 
Obviously the scattering term has only negligible effects in both limits M — > and M — > 1, where we will get the 
correlation functions of free fermions £3 ■ .. ■ 

The low-energy physics of TL a can be easily solved by the Abelian bosonizationOeS We first linearize the dispersion 
around k — ifcp (Fig. 0). We then bosonize the fermions in the upper band: 



ip(x) 



1 



-i-nMx / ao — i<*/w[(l> (x)-\-9(x)] 



(22) 



where the bosonic fields <fi(x) and 9(x) obey 
as 



(x), 6(y)\ = iO(y — x). Using these fields, we write the Hamiltonian Ti a 



dx 



g \dx 
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(23) 



where v is the Fermi velocity, and g is a parameter determined by the interaction: g < 1 (g > 1) when J± > 
(J± < 0), and g — ► 1 as M — > 0, 1. Incidentally, g is related to the compactification radius R of the field (j> by 
g = l/(47r i? 2 ). We now need to express 4>+ an d m terms of <f> and 8. Once this is done, it is straightforward to 
calculate the correlation functions since Tl a is a free-boson Hamiltonian. First we note that for states near the Fermi 
surface we have 



c \i,k c R,P + c L,fe Ci .P ~ a k a P + ^fc^P' 



(24) 



where we used the approximation (fk 
1 



Lp p w (pf. F - Using Eqs. ( |22| ) and ( |23 ) and discarding the d k fermions, we find 



- ^ + ipUx)ip L (x):fa ^(x)tp(x) = — + + — cos[2ttMx + V4tt^(x)]. (25) 

7r dx ' " ao y 7r ax 7rag 



G 



It follows that 



+ 7 1 ,„ cos(27rMa;)(cos[\/47r0(a;, t)] cosh/i^O, 0)l) a , 
where ( ) a represents the average taken in the ground state of Ti. a - The averages are found to be 

(d x <j>(x,r)d x <j>(o,o)) a =-f ( 1 + - — l^v 

47T \ (a; + i?;t) z (^ — ivt) z J 



cos[v47t0(x,t)] cos[V47T0(O, 0)]) a 



1 



2 \^a; 2 + v 2 t 2 

where ao is a short-distance cutoff of order ao . Integration of Eq. (^5|) then yields 



«0 



(x,r). 



We have neglected the contribution from the oscillating term. We thus get 

(cos[v / 7T</ , + (a:, t)] cos[^/^^(j) + (0, 0)])+ = cos(7rMx/ao)(cos[v / 7i : </ , (x, t)] cos[v / 7i : 0(O, 0)]) a 



1 / a 2 x9/4 

2 cos(7rM " /ao H^fe 



We next consider e l V^s + (x) _ p rorn gqg (17a) and (17t) we can express it as 



-iv^F0 + (a;) _ / 7ra 



(26) 

(27) 
(28) 

(29) 



(30) 



(31) 



Using the same approximation as in the derivation of Eqs. ( p4| ) and (PH), we get ipn( x ) ~ 4>( X )/V2 and i^l(x) ps 
iip(x)/y/2. Here we have made a further approximation y>fc F ps 7r/2 which is valid for < M <C 1. From Eqs. (p^), 
@, and @ we find 



*V?FM*) = |_L + sinpTrMx/ao + n/I^x) + (vr/4)] | e -^(*)+W4. 



(32) 



From Eq. (J32|) we may write the correlation functions in Eqs. (14c) and (14d) as 



e i^B+(x,r) e -i^e+(o,o) 



__ e 2irieMx/ao ( e iy/we(x,r) e ieV4^4>(x,r) e ~ieVin^(0,0) & —i v -"10.(11 



e=±l 



_ _ _ _ v —iixMx I a , _ _ _ _ . 

e i^9+(x.T) e -t^<k + (x,r) e iv^<A+(0,0) e -iv^F9+ (0,0)^ _ f ^ e i v ^i : e(x,r) e -i v ^F0(x,r) e i v ^F0(O,O) e -iv^Fe(0,0) ^ 



and 



a i-KMx/aa 



e i^d(x,r) iy/¥<t>(x,r) -iV¥0(O,O) e ~i^8(0,0) 



-37riMx /ao 



e i v ^6(^,'r) e -jv / 9iT : 0(a::T) e iv / 9iT r 0(O,O) e -iv / 7i : e(O,O) 



E 

£=±1 



,i^e + (x,T) e ieV^$+(x,T) e ieV*<t>+(0,0) e -iV^6+(a,0)^ 



e=±l 



ie7TMx/a I iy/7rO(x : T) ieyfir<f>(x ,t) — iey/Tr<f)(0,0) — i^/w$(0,0) 



The averages in the above equations are given by 



-i-Kn/2 



0() 



x + ivr 



(l/V9-«V9) 2 /4 



do 



X — IVT 



(l/Vs+«V?) 2 /4 



Combining these results together, we finally obtain the dynamical spin-spin correlation functions: 



{SZ(x,t)S5(0,0)) =m 2 -^ 



1 



1 



[x + ivr) 2 (x — ivr)" 2 



(S*(x,t)S^(0,0)) = C 2 (-l) x cos(TTMx)K (mr/v) 

i/4g 



+ Ci cos(27rMa;) 
s/4 



o 



X 2 + V 2 T 2 



'0 



(5+(x J r)5-(0,0)) = C- \Y{C, 



.T 2 + V 2 T 2 



X 2 + V 2 T 2 

c 4 



,-,2 



87T 2 f V ^ 2 + "5 2 T 2 



±- — —Koymr /v) 



X 2 + V 2 T 2 

2d e ±l7:Mx 



(i/2V?-V9) 2 



-2 _-2iriMx 

a e 



a e 



^i-xMx 



Att 2 v 



r,2 



T 2 _|_ {}2 T 2 



(l/y9-V9) 2 /4 



a; + tot 



(x + tot) 2 (x — ivr) 2 



d-Ko(mr/v) + (x — ► —a;) 



(33) 

(34) 
(35) 
(36) 



«o 



a e 



a; + tot 



a; — ivr 



(37) 



where C's are positive numerical constants, and we have set clq = 1. The correlation function (S~(x, r)$+(0, 0)) is 
obtained by replacing M with — M in Eq. (^). In Eq. ( |37| ) we have discarded the term proportional to e ±37TlMx 
decaying much faster than the kept terms. We should therefore regard Eqs. (|34| ) (|37|) as listing only the leading 
terms. As noted in Sec. II, we may expect that (S„(x, t)S^(0, 0)) should contain algebraically decaying terms that 
are proportional to cos(27r/Ma;) with any integer I. From Ref. [lg, we expect that the correlator (Sq(x)Sq(O)) should 
also have terms proportional to cos(27r/Ma;) which decay as x~ 21 9 . The appearance of the term proportional to 
e -3mMx SU gg es t s that (Sq(x,t)S^(0 : 0)) should have terms proportional to K (mr/v) cosjn(2l + l)Mx]. We note 
that the equal-time correlations (S§ (x, 0)S§ (0, 0)) and (S^ 1 (x, 0)5J(0, 0)) agree with Eqs. (|7a| ) and ( |7b| ) if we identify 
rj with 2g. As is well known, the strongest correlation is (S^(x, 0)SJ(0, 0)) ~ (— l) 3 ^ -1 /' 7 . Note also that the value of 
the exponent is consistent between the weak- and strong-coupling approach: g < 1 (g > 1) for J± > ( J± < 0) and 
g — > 1 as M — > 0, 1. Another interesting finding is that the exponentially decaying equal-time correlation functions 
have different phases by tt/2: (5£(x,0)S£(0,0)) cx cos(TrAfx) and (S$ (x, 0)S<T(0, 0)) cx sin(TrMaj). 
Now we are in a position to calculate the dynamical spin structure factors defined by 

-i /*oo />oo 

S ? f = — / dx dt(Sl(x,r = it + 0+)5^(0,0))e-^+- t , (38a) 

where t is a real time and the correlation functions in the real time are obtained by replacing r — > it + + . We may 
also calculate it from 



= -Im .ton / <fc / dr(5«(a:,r)5^(0,0)) e -^ 

for w > 0. From the obvious relation Sg"(q,(j) — Sg°(—q,u>) = S""(q + 2tt,uj), we assume < q < tt 
following discussion. We first consider the correlation functions showing the quasi-long-range order. Using Eq. (Al) 
in Appendix, we get from Eqs. (p4j) and (J36|) 



(38b) 



the 



S Z » Z {qM 



2TrM 2 S(q)d(uj) + -£^6(wWw - vq) 
2ttv 

C ~ 2 



4v 2 /5 2 



7rC 3 a5 



5[r(l/4 5 )] 2 
TrC^a 2 , 



0(w + v(q — tt))Q(lo — v(q — 7r)) 



lo 2 -i> 2 {q-2irM) 2 

,~,2 1x2 \ l-l/4ff 



1-5 



(39) 



4v 2 /d 2 



lo 2 — v 2 {q — tt) _ 
6(w - w[ g - tt(1 - 2M)])0(w + v[q - tt(1 - 2M)]) 
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u> - v[q - 7r(l - 2M)] 
2v/do 

-6(w - v[q - tt(1 + 2M)])9(< 
oj + v[q-ir(l + 2M)]\ 



2v/da 



2v/do 



oj + v[q-n(l -2M)) 
v[q + n{l+2M)]) 

2v/do 
- v[q - n{l + 2M)} 



J7+-1 



i-n- 



1-17- 



(40) 



where 7/± = [(l/2 v /g) ± y/g] 2 . The structure factor S , ~ + (q, w) can be obtained from S+~(q,ui) by M — » — M. Being 
obtained from the long-distance asymptotic expansions ([34|) and (p6|), each term in Eqs. ([39|) and ( pp[ ) describe the 
behavior of the structure factors correctly only near its low-energy threshold. For example, the last term (oc C\) in 



Eq. ( p9| ) is valid only for uj -C nvM 
structure factors are shown in Figs. 



q — 2nM\ <C nM) and cannot be extended to \q\ < nM. The supports-pf these 
They are essentially the same as those of the S = | XXZ chainOJ except 
that the small M in the ladder corresponds to the nearly polarized state in the S = A chain through the relation 

M = 5 + (S' 2 )- The strongest divergence is at q — n of (g, a;) : S^(q,u>) oc [cj — v(q — n)]~ 1+1 / 49 . The exponent 
approaches —3/4 as M — > 0,1. We note that the boundaries of the supports of these structure factors became all 
straight lines because of our linearization of the dispersion relation in the continuum limit. This is an artifact of the 
approximation, and the true boundary lines should be given by some nonlinear functions. Furthermore, some of the 
boundary lines may be parts of a single curve. 

We next consider the massive components. Using Eqs. ([sB]), ( |38b| ), and (A3), we get 



S**(q,uj) = C 2 vv 



«o 



9/2 



2v) [T{g/A)f 
(M->-M), 



,1k 



e(w - e(q - k - 7T + 7rM) - v\k\) 



e{q-k-n + nM){[uj - e(q - k - n + nM)} 2 - C 2 fc 2 }i-9/ 4 



(41) 



2„2 



The minimum energy above which S zz (q,uj) > is ui = e(q — 7r(l — M)) around 



where e(q) = \po"q 

q = 7r(l — M). Near this threshold energy the structure factor reduces to 



S* z (q,Uj)=C 2 



TTdQV 



2mvT(g/2) 



Q(u - s(q - n + nM)) 



v/d 



uj — e(q — n + nM) 



l-ff/2 



+ (M^-M), 



(42) 



where v\q— n(l =F M)\ <C mv/v and \ui — e{q — n(l^pM))\ <C m(v/v) 2 are assumed. The support of S zz (q,uj) is shown 
in Fig. ||. We see that S zz (q, w) diverges at the low-energy threshold as S zz (q, uj) oc [uj — e(q — n + nM)]~ 1+9 / 2 . Since 
the two thresholds uj = e(q — n±nM) intersect at q — n, we expect to have a peak at q — n and uj — e(nM) = H . The 
exponent approaches —1/2 as M — > 0, 1, and the singularity is even stronger for < M < 1, where g < 1. Note that 
the exponent jumps from —1 to —1/2 when H crosses H cl from below. In the strong-coupling limit of the ladder, the 
square-root divergence may be understood in the following way. The correlation function (S z i S z j) is a propagator of 
the ti t o bosons. If we ignored the interaction with the ti_ + bosons, we would get the massive free-particle propagator, 
K Q (mr /v). Due to the interaction the motion of the t i boson is necessarily accompanied by a superfluid flow of the 
ti t + bosons. Its main effect in the low-density limit amounts both to multiplying the free propagator by that of the 
hard-core 7^ + bosons oc (x 2 + v 2 t 2 )~ 1 / 4 and to shifting the momentum by nM. The Fourier transform of the product 
has the square-root divergence at the threshold. 

Finally we consider (q, uj). This can be obtained from Eqs. (p7| ) and ( |38b| ) as described in Appendix. The 
result is 



S^ T (q^) 



«o 



AnT( m )T{ m + 1) \2vJ 
+ (?-»-«), 



2)/() 



dk 



Q(uj - e(q - k T nM) - v\k\) [uj - e(q - fc =F ttM) - vk] 



'70 



e(q - k =F nM)[uj -e(q-k=f nM) + vk} 1 - 
3 1 

-e(q - k =F nM) - -v(q - k =p nM) ± m 



no 



where i]o = [(l/^/g) — v / g] 2 /4. Near the lower edge < uj — e(q =p nM) -C m, it may be approximated by 

5a 



s . r J^°(» { U - £{q -* M) ) 2V ° + {M^-M). 
8?r7;r(277o + 1) V u / a o / 



(43) 



(44) 
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In general the exponent 2rjo > and approaches as M — > 0, 1. The support of SQ~(q,uj) is shown in Fig. [| Like 
S* z (q,uj), Sq (q,v) has a peak at (q, w) — (0, if), where the two thresholds w = e(q± ttM) cross. The structure 
factor SQ + (q, uj) is approximately equal to Sq~ (q,uj)/5 near the low-energy threshold, u > sfq^fTrMj) and q w ±7rM. 
That SQ + (q,ui) is much smaller than Sq (<?, w) is consistent with the result in the strong-coupling limit where the 
ti t o bosons are absent in the ground state [see Eqs. @ and @]. 



IV. DISCUSSION 



We shall discuss implications of the results we obtained for the Heisenberg ladder to the 7r-resonance mode in the 
SO (5) symmetric ladder model. As pointed out in Ref. 0, there is an analogy between the quantum phase transition 
driven by the chemical potential in the SO (5) symmetric ladder model and the field- induced phase transition in the 
Heisenberg ladder. Obviously, thcpehcmical potential plays the role of the magnetic field. The analogy is most clearly 
seen in the strong-coupling limit .113 At half filling the ground state of the Eq phase or the Mott insulating phase 
discussed in Ref. [17| is a state in which all the rungs are in the spin singlet state. When the chemical potential 
is zero, there are fivefold degenerate low-lying massive modes above the ground state. The five modes consist of 
a S = 1 magnon triplet, a hole pair state where two holes are placed on a single rung, and a state where two 
additional electrons are put on a rung. When the chemical potential is turned on, the energy of the hole-pair 
(electron-pair) excitation decreases (increases) while the magnon triplet is not directly affected by the chemical 
potential. Thus, we see that the hole pair corresponds to the S z = 1 magnon or the t it+ boson in the Heisenberg 
ladder. The triplet magnon in the SO (5) model is an analog of the S z = magnon (t^o boson) in the Heisenberg 
ladder. Furthermore, the low-energy effective Hamiltonian for the hole-pair excitations in the strong-coupling limit 
is similar to the effective Hamiltonian for the bosons. That is, hole .pairs may be viewed as hard-core bosons 
which repel each other when two hole pairs sit on neighboring rungs .E3EIJ Let us find operators playing the role 
of the spin operators in the Heisenberg ladder. First, the operator corresponding to S z i should change a singlet 

rung |f2) = (1/V^)(ct*<it . — c||d||)|0) into a hole pair or |0). Here cj a and d\ a are creation operators of an 
electron with spin a on the zth rung of upper (c) and lower (d) chains. Obviously the ci-wave pair operator Aj is 
such an operator: A, = (c^jtiij — Ci^di^)/\/2. Second, the operator corresponding to S§ i ~ t\ + ti, + should be 
a number operator of hole pairs. It is given by N{ = 1 — (1/2) E ff ( c !ij c i ! '' + ^Jct^.c)- Finally, from the relation 
(S^ c — S^~ d )\Q,) = (c\fCi t i — d\^di ; i)\Q.) = — \/2cJ|d| j|0), we find that S^ c — Si t d creates a triplet magnon from a 
rung singlet. Thus we conclude that the spin operator Si_ c — Si_d is an analog of S z { . 

When the chemical potential is increased beyond the charge gap which equals the spin gap in the SO (5) symmetric 
model, the ladder is doped with the charge carrier (holes) and becomes superconducting with the d- wave-like symmetry. 
For spatial dimension greater than or equal to two, the superconducting order is long-ranged in the ground state, and 
this gives rise to a (5-function peak or the 7r-resonance in the dynamic spin structure factor. In one dimension , . however , 
the order is quasi-long-ranged, and therefore the peak is expected to be replaced by a power-law singularityOGj The 
threshold energy at q — n is also shown to be equal to the chemical potential.lHl These features are readily reproduced 
from our results for the Heisenberg ladder model. 

From the approximate mapping we discussed above, the ci-wave pair correlation function is expected to show the 



quasi-long-range order corresponding to the XY order in the Heisenberg ladder, Eq. (7b) 



where the exponent rj is presumably smaller than 2 and approaches 2 in the limit where the hole density S goes to 
zero. The correlator has no (— l) 4 "- 7 factor because the hole-pair mode has a minimum energy at q = 0. The charge 
density correlation is related to {Sq^Sqj}, Eq. ([fa]) , and is also quasi-long ranged: 

(N iNj ) - 6* oc COS( , 27r ^- j|) = C0S( f^~ j|) , (46) 

where we have used the relation between the hole density and the Fermi wave number, 7r(l — S) = 2kp. The 
result that the correlations of the d-wave superconductivity and AUf charge density wave show pojwer-law decay with 
the exponents whose product is 1 was also obtained by Nagaosa for a generic two-chain model.c3 Finally the spin 
correlation of the SO (5) ladder is expected to be 

" S&X^c " S ld)) <* cosM|z - j\) K f l _~^\ (47) 
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where £ is the correlation length determined by the spin gap. The spin structure factor is then 

S{q,u) oc {uj-e(q-n±nS)}- 1+i >/\ (48) 

where e(q) is the magnon dispersion at S = 0. From Fig. [5] we see that the threshold energy at q — ir is determined by 
the chemical potential, as expected. Although the exponent 77 depends on the detail of the model, we can generally 
conclude that it is 2 in the low-density limit of holes (S — > 0), where we may regard the. hard-core bosops as free 
fermions (g = 1). This universal exponent was independently found by Ivanov and LeeEj and by SchulzLj for the 
t-J ladder and was also obtained by Konik et al. for the SO(8) Gross-Neveu model. Ea We thus find that the spin 
structure factor has a universal square-root divergence at the critical point. When the superfluid density is finite, 
the interaction between bosons becomes important and modify the exponent, as we saw in the Heisenberg ladder 
model. We expect that the square-root singularity is a universal feature for the spectral weight of a gapped excitation 
generated by injecting a massive particle to a superfluid in the low-density limit. 

Although our argument above is based on the analogy and approximate mapping, our results should be valid as 
long as the weak-coupling and strong-coupling limits are in the same phase. We notice that Eq. ( ff7f ) is the same 
as the "mean-field" result given in Sec. VII of Ref. ^0|. The validity of this result is, however, questioned by Lin et 
al. as it misses the existence of the bound states such as the Cooper pair-magnon bound states found in the SO (8) 
Gross-Neveu model. On the other hand, we didn't find such a bound state in our weak-coupling calculation. It is not 
clear at the moment whether this is due to the approximation we made, for example, concerning the interaction term 
in Tt± . It was shown by Damle and SachdevcH that this term can indeed lead to a bound state of two magnons when 
H = 0. The fate of the bound state in the gapless phase is an open question. 

Finally we conclude this paper by summarizing our results on the spin correlations in the gapless phase of the 
two-leg Heisenberg ladder. We have obtained the dynamical spin-spin correlation functions and the structure factors, 
extending the bosonization theory of Shelton et al. to the gapless regime. The correlation functions are classified into 
two categories: algebraically decaying ones, (S§ (x, t)Sq (0, 0)) and {S^(x, t)SJ(0, 0)), and exponentially decaying 
ones, (S*(x, t)S*(0, 0)) and (S^(x, r)S' ( ^(0, 0)). We have also found that the terms oc cos(2irlAlx) (i. integer) are 
quasi-long-ranged, while the terms oc cos[(2^ + l)irAlx] are short-ranged. The exponents of the correlation functions 
are controlled by the single parameter g, which is smaller (larger) than 1 for J± > (J± < 0). The parameter g 
approaches 1 in the limits M — > 0, 1. The structure factors have power-law singularities at the lower edges, and the 
strongest divergence is at uj — ±v(q — ir) in S^(q, ui) due to the dominant XY spin correlation. The next strongest 
singularity is found at the lower edge of S* z (q, uj): uj = {v 2 [q — 7r(l ± M)] 2 + m 2 } 1 / 2 . The exponent is universally 
given by —1/2 in the limits M — > 0, 1. 
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APPENDIX: INTEGRALS 



In this appendix we list integral formulas we used to calculate the dynamical structure factors. 
For the gapless modes, we need the following integral: 

„ — iqx-\-iut 

dx I dt- 



(x +vt- i0+)T+ (x -vt + i0+)T- 



= 9(c-^)9(c + ^) 2 ! e (-*L-) (_JL.) . (Al) 

For the structure factor S* z (q,u) we first take the Fourier transform of the correlation function in the imaginary 
time: 



1 , (X 2 +V 2 T 2 )1 



DO 



(25)27-1 r( 7 ) (uj-is) 2 +e 2 (q-k) 
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v r(i-7) 
(2i) 2 ^ 1 r( 7 ) 



dk 



e{q-k){[uj + ie{q-k)} 2 +v 2 k 2 } 1 -f J m (Q - iv) 2 + e 2 {q - k) 



00 ^2sin(7r 7 )(i/ 2 ~i 2 k 2 y- x 



(A2) 



After the analytic continuation we take the imaginary part to find 

1 ' &(oj-e(q-k)-v\k[) 



ImI* z {q,u + iQ + ) 



(2«) 2 7-i vr(7) 



dk- 



(q - k){[w - e{q - k)} 2 - C 2 fc 2 }!-7 



(A3) 



for uj > 0. When v\q\ <C mv/v, e(q — k) + v\k\ w e(g) + {}|A;|. In this case we may approximate the last integral as 



m 



e(o; - e(g) - g|fc|) , . gQyjyg) 27 _ x 

{[w - e(g)] 2 - v 2 *; 2 } 1 ^ v K " vm L v " 



where B(a,b) is the beta function. From Eqs. (|A3|) and (A4) we finally obtain 



Iml* z {q,uj + i0+) = Q(u-e{q)) 



r(27) mv \uj — e(q) 



1-27 



(A4) 



(A5) 



which is valid for uj — e(q) <C m. 

We next consider S^(q,ui). According to Eq. (p^), we need the following Fourier transform: 



— 00 o — 00 



' ' (x z + v 2 r z y (x — IVT) 

° dk [°° d i v ^ ^k^ 1 (y — ivk) 1 
(uj — v) 2 + e 2 (q — k) 



.r(i-7) v 



r(l+ 7 ) (2v) 2 -r J_ 

r(i- 7 ) v 



00 •/ —00 



dk 



Z7r[a; + i5fc + ie{q — k)} 1 



r(l+ 7 ) (2w) 2 t J_ 00 \s(q - k)[w - ivk + ie(q - k)} 1 -^ 
After the analytic continuation we obtain 



2 sin(7r 7 ) / dv 



(y + vk) 1 {v — vk) 



7-1 



5|fe| {uj - iv) 2 + e 2 (q- k) 



lml+-{q,uj + i0 + ) = 



7(2^ Vr(7) 



dkQ(u - e(q - k) - v\k\) 



[uj — e(q — k) — vk] 1 
e(q -k)[u- e(q - k) + vk} 1 —' ' 



Using the same approximation as in Eq. (A4), we obtain 

imtf-( ffl w + *o+) = e( w - e(ff) ) ^^pr^ ( ^ ) 

v 7 1 (1 + 27) tow \ v ) 

for < uj — s(q) <C to and <C mv/v. In the same way we get 



27 



(A6) 



(A7) 



(A8) 



Ihi/q ^{q,uj + i0 + ) = Im lim z / / dr 



OO /"CO 



U / 7T 



00 «/ —00 
2 /» 00 



Ko(my/x 2 + v 2 r 2 /v) 
(x 2 + £ 2 r 2 ) 7 (a; + zC>t) 



dfce(w-e(g-fc)-u|fe|) 



7 (2«) 2 -/ Vr( 7 ), 

which reduces to Eq. ( |A8| ) for < uj — e(q) <C to and -C mv/v. 



[uj — e(q — k) + vk] 7 
e{q-k)[uj-e{q-k)~vk} 1 -<'' 



(A9) 
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FIG. 1. Schematic picture of the upper and lower bands: E(k) — ±v v 2 k 2 + m 2 — H. The negative-energy states are filled. 



The long-distance behavior of the correlation functions are determined by the low-energy excitations around \k\ = 
the dispersion is linearized. 
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2nM q 

FIG. 2. Support of So z (q,ui). The shaded region shows where S(j z (q,uj) is nonzero. 
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FIG. 3. Support of St (q,<+>)- The shaded regions show where 5^ 
uj — ±v(q — n). The next strongest singularity is at uj = ±v[q — n(l ± 2M)]. 



is nonzero. The strongest divergence is at 
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FIG. 4. Support of S 7r + (q,ui). The shaded regions show where S 7r + (q, uj) is nonzero. The strongest divergence is at 
uj — ±v(q — 7r). The next strongest singularity is at uj = ^v[q — n(l ± 2M)}. Note the difference from Fig. [| 
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FIG. 5. Support of S* z (q,io). The shaded regions show where S zz (q,u>) is nonzero. 
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FIG. 6. Support of S^~{q,uj) and S'+^u;). The shaded regions show where S^(q,uj) are nonzero. 
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